home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / clags2.z / clags2
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCLLLLAAAAGGGGSSSS2222((((3333FFFF))))                                                          CCCCLLLLAAAAGGGGSSSS2222((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CLAGS2 - compute 2-by-2 unitary matrices U, V and Q, such that if ( UPPER
  10.      ) then   U'*A*Q = U'*( A1 A2 )*Q = ( x 0 )  ( 0 A3 ) ( x x ) and  V'*B*Q
  11.      = V'*( B1 B2 )*Q = ( x 0 )  ( 0 B3 ) ( x x )  or if ( .NOT.UPPER ) then
  12.      U'*A*Q = U'*( A1 0 )*Q = ( x x )  ( A2 A3 ) ( 0 x ) and  V'*B*Q = V'*( B1
  13.      0 )*Q = ( x x )  ( B2 B3 ) ( 0 x ) where   U = ( CSU SNU ), V = ( CSV SNV
  14.      ),
  15.  
  16. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  17.      SUBROUTINE CLAGS2( UPPER, A1, A2, A3, B1, B2, B3, CSU, SNU, CSV, SNV,
  18.                         CSQ, SNQ )
  19.  
  20.          LOGICAL        UPPER
  21.  
  22.          REAL           A1, A3, B1, B3, CSQ, CSU, CSV
  23.  
  24.          COMPLEX        A2, B2, SNQ, SNU, SNV
  25.  
  26. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  27.      CLAGS2 computes 2-by-2 unitary matrices U, V and Q, such that if ( UPPER
  28.      ) then
  29.            ( -CONJG(SNU)  CSU )      ( -CONJG(SNV) CSV )
  30.  
  31.        Q = (     CSQ      SNQ )
  32.            ( -CONJG(SNQ)  CSQ )
  33.  
  34.      Z' denotes the conjugate transpose of Z.
  35.  
  36.      The rows of the transformed A and B are parallel. Moreover, if the input
  37.      2-by-2 matrix A is not zero, then the transformed (1,1) entry of A is not
  38.      zero. If the input matrices A and B are both not zero, then the
  39.      transformed (2,2) element of B is not zero, except when the first rows of
  40.      input A and B are parallel and the second rows are zero.
  41.  
  42.  
  43. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  44.      UPPER   (input) LOGICAL
  45.              = .TRUE.: the input matrices A and B are upper triangular.
  46.              = .FALSE.: the input matrices A and B are lower triangular.
  47.  
  48.      A1      (input) REAL
  49.              A2      (input) COMPLEX A3      (input) REAL On entry, A1, A2 and
  50.              A3 are elements of the input 2-by-2 upper (lower) triangular
  51.              matrix A.
  52.  
  53.      B1      (input) REAL
  54.              B2      (input) COMPLEX B3      (input) REAL On entry, B1, B2 and
  55.              B3 are elements of the input 2-by-2 upper (lower) triangular
  56.              matrix B.
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCLLLLAAAAGGGGSSSS2222((((3333FFFF))))                                                          CCCCLLLLAAAAGGGGSSSS2222((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      CSU     (output) REAL
  75.              SNU     (output) COMPLEX The desired unitary matrix U.
  76.  
  77.      CSV     (output) REAL
  78.              SNV     (output) COMPLEX The desired unitary matrix V.
  79.  
  80.      CSQ     (output) REAL
  81.              SNQ     (output) COMPLEX The desired unitary matrix Q.
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.